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Abstract 


This paper deals with a novel proof of convergence of He’s variational iteration 
method applied to nonlinear partial differential equations by proposing a new for- 
mulation for this technique. 
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1 Introduction 


Recently He [8] has written a survey article and some new asymptotic tech- 
niques with numerous examples. The limitations of traditional perturbation 
procedures are illustrated. Various modified perturbation techniques are in- 
troduced, and some mathematical tools such as variational theory, homotopy 
technique, and iteration technique are proposed to overcome the shortcom- 
ings. For the nonlinear oscillators, all the reviewed schemes produce high 
approximate periods, but the accuracy of the amplitudes cannot be amelio- 
rated by iteration. The emphasis of this author [8] is on the variational ap- 
proaches, parameter-expanding methods, parameterized perturbation tech- 
nique, homotopy perturbation method, iteration perturbation procedure and 
ancient Chinese methods. Variational approaches to soliton solution, bifur- 
cation, limit cycle, and period solutions of nonlinear equations including the 
Ritz method, energy technique, variational iteration method are illustrated 
in his paper [8]. 

The variational iteration method (VIM) plays an important role in recent 
researches in this field. This method is proposed by He [4, 6, 5, 7] as a 
modification of a general Lagrange multiplier method [11]. It has been shown 
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that this procedure is a powerful tool for solving various kinds of problems 
(e.g., see [1, 3, 13, 12]). 

In this work, we adapt the technique to nonlinear partial differential equa- 
tions and we prove the convergence of this method by proposing a new for- 
mulation of the method. 


2 The Variational Iteration Method 


The idea of VIM is very simple and straightforward. To explain the basic idea 
of VIM, we consider an one dimension general first order nonlinear partial 
differential equation as follows with the assumption that the equation has the 
unique solution (note that one can consider a nonlinear partial differential 
equation in higher dimensions and also more general form without loss of 
generality): 

P(t, v7, U, Uy, Ue, ***) = 9, (1) 


with the specified initial condition (i.e., u(x, to) = f(x)). We assume that the 
nonlinear operator ® is continuous with respect to its arguments and u(z,t) 
is an unknown. We first consider Eq. (1) as follows: 


A[u(a, t)] + N[u(a, t)] = 0, (2) 
with, for example, the assumption a(x,t) 4 0 
Alu] = a(a, thy + B(a, thu 
and 
N[u] = ®(t, x, u, uw, ue, ++) — a(x, t)uy — B(a, t)u, (3) 


where, as shown above, A with the property Ay = 0 when y = 0 denotes the 
linear operator with respect to u and N is a nonlinear operator with respect 
to u. We then construct a correction functional for Eq. (2) as [7]: 


Unqi(@,t) = Un(a,t) +f Aww,t,s) { O(2, 8)tin, (£2, 8) + B(x, 5)tn (2, 8) 
+Nitin (x, s)] } ds, (4) 


where uo(a,t) is the initial guess and the subscript n denotes the n-th itera- 
tion, and A(z,1,s) 7 0 denotes the Lagrange multiplier, which can be identified 
efficiently via the variational theory, and u,, is considered as a restricted vari- 
ation [7], ie., du, = 0. 

Taking variation with respect to the independent variable u,,, noticing that 
dUn(x, to) = 0 and by making the correction functional stationary, we obtain 
bUn+1(@,t) = 0 and therefore we have the following: 
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t 


bUn4i1(@,t) = dUn(a,t) + 5 | Aa,t,s) { O(2, 8)Un, (2, 8) + B(x, 8)tn(z, 8) 
to 
+Nitin (x, s)]} ds 
= bUn(2, t) T a(x, 8)X(x,t,8)OUn (%, s) | at 


-f 7 (a(x, 8)Xa,t,8)) — Be.) Xcoa } Un (x, )ds 
(il oe aan) 6un(Z, 8)|,—4 

_ [ ie (cule, 8) Avea,e) — A(t, s)Aena} Sun (a, 8) ds 
= (5) 


Therefore, we have the following stationary conditions: 


(2, 8) Nts 2.2) | py = —-1, (6) 


0 
3s (a(z, 8) (x,t,s)) — B(@, 8)A(a,t,s) = 9. (7) 


Hence, the Lagrange multiplier can be readily identified as 


Ne,t,) = ae exp ([ eh) 79) ar) (8) 


a(a,T) 


As a result, we have the following variational iteration formula: 


t 

Un4i(@,t) = Un (a, t) + A (ait a) OS, ©, Ug BS), Uns (a, 8), Uy, (By S) 5° ** ) ds. 
to 

(9) 


Accordingly, the successive approximations u,(a,t),n > 0 of VIM will be 
readily obtained by choosing all the above-mentioned parameters. Conse- 
quently, the exact solution may be obtained by using 


u(a,t) = Jim Un(a, t). (10) 


The initial guess can be freely chosen with possible unknown constants, 
it can also be solved from its corresponding linear homogeneous equation 
A[uo(#,t)] = 0. It is important to note that for linear problems, the exact 
solution can be obtained easily by only one iteration due to the fact that the 
auxiliary function can be suitably identified [9]. For nonlinear problems, in 
general, one iteration leads to highly accurate solution by VIM if the initial 
solution is carefully chosen with some unknown parameters. 
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3 Convergence Theorem 


The variational iteration formula makes a recurrence sequence {u,,(x, t)}. Ob- 
viously, the limit of the sequence will be the solution of Eq. (1) if the sequence 
is convergent. In this section, we give a new proof of convergence of VIM in 
details by introducing a new iterative formulation of this procedure. Here, we 
suppose that for every positive integer n, Un € Clto, T], and {un, (az, t)} and 
{un, (x, t)} are uniformly convergent. 


Lemma 3.1. The variational iteration formula (9) is equivalent to the fol- 
lowing iterative relation 


Alun+i (2, t) — Un(a, t)] = —®(t, x, un (ax, t), Un, (@, t), Un, (@,t),--+), (1) 


where A is as noted in (3). 


Proof. Suppose u,, and u,+1 satisfies the variational iteration formula (9). 
Applying Ba into both sides of (9), results in 


0 : On x,t,s O 
Sjltins (#4) — un (2, )] =| 2) bds+ S[Netalace (12) 
0 
Now, by using conditions (6) and (7), and Pgs) = nea we will have 


a(a, t) 5; [Un+1 (x,t) — Un(2, t)] + B(x, t)[un41(2,t) — Un(2, t)] = 


Te a eo 


From the definition (3) of A, we obtain 


Alun+i(a,t) — Un(ax, t)] = —®(t, x, un (x,t), Un, (@,t), Un, (x, 6),---). (14) 


Conversely, suppose up and Un+1 satisfies (11). Multiplying (11) by \(a2,4,s), 
in view of the definition of A and \(24,3) # 0, and next by integrating on 
both sides of the resulted term from to to t, yields 


ie A(a,t,8) (2; [gener (® s) _ Bun (a, s)] ds + i A(a,t,s) B(@, 8) [Un4i (2, s) _ Un(a, s)| ds 
=— i. A(x,t,s) P(s) ds. 
(15) 
Using simple integration by parts, the expression (15) becomes 


a(x, t)Ax,r,t) lUn41 (x, t) — un(a, t)] — i, (2 (a(a, 8)Xa,t,s)) — B(®s 8) A(a,t,8) [ni (2, 8) 
—un(a, s)] ds = — iia A(wx,t,s) &(s) ds, 
(16) 
which exactly results (9) upon imposing conditions (6) and (7), ie., 


t 
Un+i(a@,t) = Un(a,t) + / A(z,t,2) 2 (8, £, Un (2, 8), tin, (2,8); Un, (©, 8),-*> )ds. 
t 
: (17) 
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and this ends the proof. 


Theorem 3.1. If the sequence (10) converges, where un(a,t) is produced by 
the variational iteration formulation of (9), then it is the exact solution of 
the equation (1). 


Proof. If the sequence {u,,(x,t)} converges, we can write 


u(x,t) = dim Un(a,t), (18) 
and it holds 
v(a,t) = lim Un+1(2,t). (19) 


Using expressions (18) and (19), and by the definition of A in (3), we can 
easily gain 


lim A [Un4i(a, t) — Un(ax, t)] = A lim [un+i(@, t)—Un(a#,t)}=0. (20) 


noo 
From (20) and according to the Lemma 3.1, we obtain 
A lim [un+i(#,t) — un(a,t)] = — lim ®(t,2,un(2,t), Unt (x,t), Ung (x, t),---) =0, (21) 
which gives us 


lim ® (t, 2, un(,t), Un, (x, t), Un, (@,t),+-+) = 0. (22) 


noo 


From (22) and the continuity of ® operator, it holds 


lim ®(t, 2, Un (x,t), Un, (7, t), Un, (@,t),+-+) 
nN—-Ooo 


= P(t, w, jim, Un(z, t), jim, Unt (x, t), im, Une (x, t), 7s ) 
= (4. fim aC.) (tat) » (im ast) oo) 
= O (t,x, u(a, t), ve(x, t), ve(z, t), +++). (23) 


Now, from Equations (22) and (23), we have 
® (t,2,0,U4,Ur,°°:) =0, to <t<T. (24) 


On the other hand, using the specified initial conditions and the definition of 
the initial guess, we have 


v(x, to) = lim Un(@,to) = f(x), since u(x, to) = Un(x,to) = f(x), n>0, 

(25) 
Therefore, according to (24)-(25), v(#,t) must be the exact solution of the 
equation (1), this ends the proof. 


Note that the above theorem is valid for the linear operator A defined by 
(3). This convergence theorem is important. It is because of this theorem we 
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can focus on ensuring that the approximation sequence converges. It is clear 
that the convergence of the sequence (10) depends upon the initial guess 
uo(a,t) and the linear operator A. Fortunately, VIM provides us with great 
freedom of choosing them. Thus, as long as uo(a,t) and A are so properly 
chosen that the sequence (10) converges in a region to < t < T, it must 
converge to the exact solution in this region. Therefore, the combination of 
the convergence theorem and the freedom of the choice of the initial guess 
uo(a,t) and the linear operator A establishes the cornerstone of the validity 
and flexibility of VIM. 


4 An Illustrative Example 


In order to illustrate the efficiency of the VIM described in this paper, we 
present one example. 
Example A much-considered example is the Burger’s equation [12,13] 


Ut + UUg — Uge = 0, (26) 


This equation was only intended as an approach to the study of turbulence 
because it exhibited some essential characteristics of the more realistic (and 
difficult) equations. This equation involves nonlinearity, dissipation, and is 
relatively simple. The VIM solves much more difficult systems. We consider 
it now to show the simplicity of a proper solution. 

According to the VIM procedure, (9), one can obtain the following varia- 
tional iteration relation (A(z,4,s) = —1): 


Unti(a,t) = un(et)— f {Un, (©, 8) +Un(Z, 8)Un, (@, 5) —Un,,, (@, 8)} ds. (27) 


The problem is completely defined when the initial condition is specified. If 
we specify u = x when t = 0, we have 


uo(a,t) = x, 

ui(a,t) = a2 — at + O(2°), 

uo(a,t) = 2 — at + at? + O(t>), 

ug(x,t) = 2 — t + at? — 2i*® + O(¢'), 
ua(a,t) = 2 — vt + xt? — ct? + t* + O(2°), 


(28) 


Thus, u(a#,t) = a/(1+t) which is the exact solution of (26). 


On convergence of He’s variational iteration method for nonlinear partial... 7 


5 Conclusion 


In this paper, we have given a proof of convergence of He’s variational iter- 
ation method by presenting a new formulation of He’s method for nonlinear 
partial differential equations. The main property of this method is in its flex- 
ibility and ability to solve nonlinear equations accurately and conveniently 
without decomposing the nonlinear terms, which makes the procedure very 
complex. This technique is a very powerful tool for solving nonlinear prob- 
lems. Furthermore, it gives an accurate and easily computable solution by 
means of a truncated series whose convergence is fast. 
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